Elastic enhancement factor as a quantum chaos probe 
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Recent development of the resonance scattering theory with a transient from the regular to chaotic 
internal dynamics inspires renewed interest to the problem of the elastic enhancement phenomenon. 
We reexamine the question what the experimentally observed value of the elastic enhancement factor 
can tell us about the character of dynamics of the intermediate system. Noting first a remarkable 
connection of this factor with the time delays variance in the case of the standard Gaussian ensembles 
we then prove the universal nature of such a relation. This reduces our problem to that of calculation 
of the Dyson's binary form factor in the whole transition region. By the example of systems with 
no time-reversal symmetry we then demonstrate that the enhancement can serve as a measure of the 
degree of internal chaos. 
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Excess of probabilities of elastic processes over inelastic 
ones is a common feature of the resonance compound 
nuclear reactions, electron transport through quantum 
dots, where it manifests itself as the weak localization 
effect, or, at last, transmission of electromagnetic waves 
through microwave cavities. This phenomenon, that is 
characterized quantitatively by the elastic enhancement 
factor i.e. the typical ratio F of elastic and inelastic cross 
sections, repeatedly attracted attention for decades [1, 
M. ITci \vK ITal . Based on the random matrix theory j2l 
universal formalism that allows of uniform treatment of 
all resonance phenomena of such a kind has been worked 
out in the seminal paper |3] . The scattering M X M matrix 
that describes the resonance processes is expressed as 



S(E) = I - z'A + - 



■A. 



(1) 



The Hamiltonian matrix H that describes dynamics of 
the originally closed intermediate system is supposed to 
belong to an ensemble of random matrices. This system 
gets excited and then decays after a while because its 
N » 1 eigenstates are connected to M open channels via 
random j^] rectangular N x M matrix A. The formula 
((TJ can be naturally interpreted 0, @] (see also 0, @] 
and references therein) within the concept of open 
systems whose internal dynamics is described by a non- 
Hermitian effective Hamiltonian 'H = H - jW where 
the connection to channels results in the anti-Hermitian 
contribution proportional to the matrix W = AA f . 
The complex eigenvalues of the effective Hamiltonian 
specify positions and widths of the resonances. 

All quantities of the physical meaning are obtained 
by averaging (...) over two independent ensembles 
of matrices H and A J^]. Specifically, the averaged 
scattering matrix (S ah (E)) = (S aa (E)) 5 ab fixes the trans- 
mission coefficients T" = 1 — KS flfl )| 2 that measure the 



part of the flow that spends essential time in the internal 
region. The scale of the energy dependence of the mean 
scattering matrix as well as the transmission coefficients 
is large and their energy variations can be neglected 
within the energy interval where the N resonance states 
of interest are situated. 

In the most interesting case of large number M » 1 of 
scattering channels (for the sake of simplicity we sup- 
pose that all these channels are statistically equivalent 
and have, therefore, identical transmission coefficients 
T) the enhancement factor F consists of two contributions 
of quite different nature. The first of them, 1 + 6\p, de- 
pends only on the symmetry class jS = 1 (preserved time- 
reversal (T) invariance) or /3 = 2 (broken T-invariance) 
of the correspondingensembles of the random effec- 
tive Hamiltonians f H la- HOIl . No enhancement originates 
from this contribution in the case of Unitary ensemble of 
complex Hermitian matrices. The second contribution 
is regulated by the ratio rj = tn/tw = MX ("openness") 
of two characteristic times. One of them, the Heisenberg 
time 



*H = = -23<TrG(E + iO)) , 



(2) 
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where G(E + z'O) = (E + z'O - H) -1 is the resolvent of the 
Hamiltonian H, characterizes the internal motion and is 
defined by its mean level spacing d. Similarly fHHl^l . 
the dwell time f w = j(Q)a^o, where (Q) = -- 2 j3(Tr^ 7 ) 
is the mean delay time, establishes the time scale of the 
open system in terms of the mean Wigner delay time 
il4Lll5ll . The dwell time is the time the incoming parti- 
cle spends in the internal region. The inverse quantity 
Tw = 1/fyv is nothing else than the well-known Weis- 
skopf width ifLsll . Only if the openness is small enough, 
rj «c 1, so that the dwell time appreciably exceeds the 
Heisenberg time, fjj <SC tw , the incoming particle has 
enough time to recognize the discreteness of the internal 
spectrum and therefore to perceive spectral fluctuations. 
Then an additional contribution appears JUlll] and the 



2 



enhancement factor takes finally the form 

Xoo 
dse-v [l - Bf (s)] 



(3) 



where Bf\s) is the Dyson's spectral binary form fac- 
tor J2] belonging to the symmetry class B. Comparing 
this expression with the delay time two-point correlation 
function 



<Q 2 > 

Bf (s)J cos (2nspe) 



(4) 



calculated for the case B = 1 in 1121] (later on, similar 
calculation has been perf ormed also in the case B = 2 fl3l0 
we arrive at the following remarkable relation between 
the enhancement factor and variance of the time delays: 



FW(tj) = l + 5 jn +r 1 varQ(T 1 ) 

-- 2 + <V - n Jo 



dsrVBf (s). 



(5) 



A new aspect of the old problem of the elastic en- 
hancement has been recently evoked in ref.|119ll (see 
also earlier semiclassical consideration in 12011 ') where 
manifestations of transition from regular to chaotic 
internal motion has been investigated in the framework 
of the resonance scattering theory. The character of this 
motion is controlled by a particle interaction parameter 
K. The dynamics of the intermediate system can there- 
fore be described by some transient matrix ensemble. 
It has been, in particular, numerically discovered that 
the elastic enhancement factor is quite sensitive to the 
strength of the interaction. This fact suggests that the 
enhancement factor can serve as an indicator of the 
degree of the internal chaoticity. In what follows we 
investigate this relation analytically. 

Generally, the short range universal fluctuations of scat- 
tering amplitudes are described by the (connected) S- 
matrix two-point auto-correlation function [3] 



C abcd (e) = <S fl0 (E + e/2)S ca *(E - e/2)> 



-nb 



cd * 



(6) 



While carrying out the ensemble averaging we, suppose 
throughout this Letter the decay amplitudes A c n to be 
uncorrelated Gaussian random quantities 



(KA [ ;; l ) A = 1 -b nm b ab . 



(7) 



This assumption is supported by the reasons of so-called 
"geometrical chaos" that have been argued in ref. l2lll . 
In the case of T-invariant systems B = 1 these amplitudes 
are real, A b n * = A b m . 

For the elastic enhancement factor the formula (|6) gives 
F = C aam (0)/C abab (0) . We start ensemble averaging with 



that over the amplitudes A keeping the internal Hamilto- 
nian H diagonal. It is convenient (though not necessary) 
to use supersymmetric integral representation. Then A- 
averaging can be fulfilled exactly whereupon the saddle 
point method can be used. In such a way we receive first 
of all 



(S ab (E)) = 5' 



ab 



l + i 



(8) 



where the function g(E) = ^(Trjp^)^ satisfies the equa- 
tion (m=M/N) 



1 N ' 



h=1 L 



E-E„ + 



1 + hym 



(9) 



The subscript P in the r.h.s. of eq. <(8j> implies averaging 
over all energy levels of the internal system with the 
joint probability distribution P({£}|?c); {E} = Ei,£2,...Ejv 
at a given value of the chaoticity parameter k. Depend- 
ing on this parameter, the distribution P({E}|k) changes 
from Poissonian distribution of fully independent levels 
(k = 0) to that of highly correlated levels what is typical 
of the Gaussian ensembles (k = 00). We assume also 
that the mean level density does not depend on x at all. 

In the limit of weak coupling to continuum y — » we are 
interested in, the approximate solution of eq. l[9| reads 



g(E) *iTr G^E + ^my) 



(10) 



We suppose below that the ratio m is also small. Then 
the mean scattering matrix reduces ll5l. l22l] to 



(S ab (E)) = d 



1 + x 



(11) 



In accordance with the conventional practice, we 
neglected the long range energy dependence of the 
mean S-matrix elements and set E = 0. (We will do the 
same in all later calculations.) The measuring the degree 
of resonance overlapping parameter x = (where d 
is the mean level spacing) should be small in the case 
of our interest so that the corresponding transmission 
coefficients equal T ~ 4x and, correspondingly, the 
openness is r\ = 4M x 



The tensor structure of the correlation function lO 

C abcd (e) = Ff(e)5 ab 5 cd + pf (e)(& K 6 bd + d w 5 ad 5 bc ) (12) 

follows from the T-invariance properties and rotational 
invariance in the channel space. The superscript B marks 
now the symmetry class of the limiting (k = 00) matrix 
ensemble. The enhancement factor reads therefore 



F^ = l + V + [Ff (e) /Ff(e)] 



e=0 



(13) 
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Now, in the leading order with respect to y the A- 
averaging results in 

Ff(N,M,y\K) = £(TrG(±my)TrG%my)) pcgm 
Ff{N,M,yW) = £(TrG(imy)G\im 7 )) p . 

Subsequent P-averaging is straightforward and leads to 
the expression 

d8e-*[l-Bf>m]. (14) 

Note that this ratio depends after all only on two param- 
eters: on the openness )] of the internal system and on 
the degree of chaoticity k of its dynamics. Finally we 
arrive at the expression 

dse^Bf (s|k) (15) 

o 

that extends the relations d3j |5]l) to the case of arbitrary 
value of the chaoticity parameter. 

The found result reduces the problem posed above to 

that of calculating the binary form factor B^\s\k) in 
the whole transient region < k < oo. The issue of 
transition from regular to chaotic dynamics has been 
attacked not once by different authors (see l23ll and 
references therein). The total solution has been found 
by now only in the case of the systems with broken 
time-reversal symmetry l24l. l25ll . The method used in 
ll25ll is the most convenient for our purpose. These 
authors has used the Brezin-Hikami's approach j26ll 
that allows of direct calculating the binary form factor 
we need. Below we restrict ourselves to the case jS = 2 
and will skip this superscript. 

The following two properties of the considered bi- 
nary form factor are obvious from the very beginning: 
B 2 (s\0) = and B 2 (s|oo) = (1 - s)0(l - s), jj. In the 
intermediate region the form factor is given by jzsll 



B 2 (s|k) = B 2 (s|oo) 

2 rl [2yVS+l] 



s+2y Vs+1 



y_„-Ks(s+2y Vs+1) 



(16) 



In fact, only the even part of the integrand contributes. 
The enhancement factor is entirely expressed via the 
function 



vj/^Ik) = n JT" ds e ->«(s + i)-st, 

~'"- s, 'E(2ks 3 / 2 ) 



= tj£° dse~ KS ^ +2Ks3 ' 



(17) 



where h(x) stands for the modified Bessel function 



and the function E(2ks j/ ) 



3/2 \ _ r -2KS 3 ' 2 h(2KS 3 l 2 ) 



decreases 



monotonously from one to zero when the argument 
2 KS 3/2 g rows n j s e asy to check that ^(/jlO) = 1, 
W(rj\co) = and ^(^lic) > in between. 



There are two equivalent representations of the enhance- 
ment factor: 



F(T]\K) = 1 + W(T]\K) + Ti-z 



'dtf 



1 Jo 



(18) 



and 



F(T/|K) = l + (l-e^)/»j 



(19) 



In particular, the first formula gives immediately F(rj\0) = 
2 when the second one reduces to the GUE result Fgue = 
1 + (1 - e~'i)/t] as 2 - \r\ + ... . More than that, it can be 
shown with the aid of eq. dl9t that the slope at the point 
/] = is universal: 



dF(rj\ K * 0) 



drj 



r,=C 



(20) 



Indeed, at any nonzero k contributions of the two last 
terms cancel each other when 77 — > with accuracy 
better than 77. 

Although one cannot obtain any exact explicit analytical 
formula, a number of approximate expressions can be 
derived from eqs. 0.8U19II . At that eq. lfl8l is useful when 
the internal chaoticity is weak whereas the second form 
is more convenient if the internal dynamics is close to 
chaotic. Behavior of the function W depends on interre- 
lation of the two competing parameters K and 77. If the 
first of them is small enough and the second one is kept 
finite we can take into account only few leading terms 
in the Bessel function power series. Corresponding con- 
tributions are expressed already in the terms of known 
transcendental functions. After that there are two possi- 
bilities: either expand these functions into power series 
with respect to the parameter x or make expansion over 
inverse powers of the parameter 77. In the first case coef- 
ficients of the K-expansion are polynomials in I/77; in the 
second case those of the l//]-expansion are polynomials 
in k. The two found in such a way expansions do not 
perfectly coincide. Nevertheless they match with certain 
accuracy that can be improved by taking into account a 
larger number of contributions. Substituting finally the 
estimated in such a way function W into eq. 1IT8I we 
arrive at 

r( \ \ o * , (6 + T7)* 2 (60 + 77(20 + 7]))^ _ 

F(t7|k) = 2 - - + + ... . (21) 

77 rf rf 

On the other hand, if the parameter of chaoticity is 
large, K » 1, calculations become appreciably simpler 
and integration in eq. i fTTt can be carried out with the 
help of the Laplace method. At that, it is convenient 
to utilize the presentation given in the second line of 
eq. il7t . Generally speaking, the exponential factor in 
the integrand has two maxima in the points Sq = and 
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sj = 1/8 (5 - 4/]/k + 3yjl- 8j]/k) .In the first of them the F(Ar|?7) 

...u„i„:„,. J 1 -ri , TI /„i..\ - n __j xu„ 2.0 r 



whole integrand equals one. Then ^(rjlx) ~ and the 
contribution of the vicinity of the point So in the enhance- 
ment factor is easily found to be 



F(/]|k) * 1 + 



7] + K (77 + K) 2 



(22) 



As to the second point si, the maximum of the ex- 
ponential reaches its largest possible value one when 
77 = 0, goes rapidly down with growing 77 and, after 
passing the inflection point s, = 9/16, disappears 
finally when 77 exceeds k/8 . But even in the most 
interesting case K » 77 > 1 contribution of the vicinity 
of the second maximum remains negligible. Indeed, 
opposite to the hight of the maximum neither its 
position si » 1, nor its width As * 0.4 noticeably 
depend on 77 . Therefore the slow varying factor can 

be estimated as H(2ks 3/2 ) » E(2ks 3/2 ) » 1 /V2~tzk 3/2 <k 1 . 

The Fig.l illustrates variations of the elastic enhance- 
ment factor depending on the increasing openness 77 and 
parameter chaoticity K. In particular, the domains of va- 
lidity of our approximations are demonstrated. At any 
given chaoticity k the factor F(t]\k) decreases with the 
universal initial slope -1/2 starting from the maximal 
value 2 up to some value 77 c (k) where this factor reaches 
minimum F m „, so that 



dF(rj\K) 



r\=r\ c (K) 



= 0. 



(23) 



The level correlations at the chosen value of the chaotic- 
ity parameter k are too weak to be resolved when 
openness 77 exceeds the value r] c (ic) and the enhance- 
ment factor returns to value 2 typical of the system 
with regular dynamics. In this way the enhancement 
factor conveys information on the degree of internal 
chaos. Finally, the required chaoticity parameter K(F m! „) 
is found as the root of the equation F (t] c (k)\k) = F m ,„ . 

Though we do not know the explicit form of the binary 

form factor B? (s|k) that describes considered transi- 
tion in the case of F-invariant systems, there is no doubt 
that, qualitatively, the behavior of the enhancement 
factor should be similar. It is worth mentioning that the 
evolution of this factor under transition from perfectly 
chaotic systems with broken F-invariance to F-invariant 
ones can easily be followed with the aid of the method 
developed in I27I1 . 

The further quite nontrivial development reported in 
l28Tl opens opportunity of deriving such a form factor 
for T-invariant systems also. Corresponding results will 
be published elsewhere. 



Summary We have considered the dependence of the 
elastic enhancement factor on the degree of chaoticity of 
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Figure 1. (color online) Enhancement factor F(j]\k) as a function 
of 77 for three fixed values (from top to bottom) k = (0.5, 5, 50). 
Black points - exact numerical results, blue lines (dashed) - 
small k approximation SS\\ . red lines (full) - large K approxi- 
mation (22) ; the straight green (gray) line shows the universal 
slope at the point 77 = 0. 



the internal part of an open resonance system. A general 
relation of this factor with the variance of time delays 
has been established for arbitrary degree of chaoticity. 
By this, the task is reduced to the search for the transient 
binary form factor. Generally, the enhancement factor 
F(?7|k) depends on both the chaoticity parameter k and 
the openness 77 the latter being the ratio 77 = tjj/tw 
of two characteristic times: the dwell time tyj and 
the Heisenberg time tu- Here the time fjj = *f is the 
time that is needed to resolve the pattern of spectral 
fluctuations in the system with Hamiltonian H when 
the time f w - 1 /Tyv is that the incoming particle spends 
in average inside this system. Only if this particle is 
trapped inside for sufficiently long time it can carry 
information on the internal chaos. Otherwise the dif- 
ference between regular and chaotic internal dynamics 
cannot be resolved. 

In this Letter, the problem posed has been thoroughly 
studied numerically and analytically in the case of sys- 
tems with no time-reversal symmetry. We showed in 



particular that the slope — ^ — - = — | remains invari- 
able for arbitrary degree of internal chaos. The recovery 
of the maximal value of F(?7|k) when the openness 77 ex- 
ceeds some value 77 c (k) that is clearly seen in the Fig.l 
is in perfect agreement with the physical argumentation 
stated in the previous paragraph. 
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